Self-ratcheting Stokes drops driven by oblique vibrations 
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We develop and analyze a minimal hydrodynamic model in the overdamped limit to understand why a drop 
climbs a smooth homogeneous incline that is harmonically vibrated at an angle different from the substrate 
normal [Brunet, Eggers and Deegan, Phys. Rev. Lett. 99, 144501 (2007)]. We find that the vibration component 
orthogonal to the substrate induces a nonlinear (anharmonic) response in the drop shape. This results in an 
asymmetric response of the drop to the parallel vibration and, in consequence, in the observed net motion. 
Beside establishing the basic mechanism, we identify scaling laws valid in a broad frequency range and a flow 
reversal at high frequencies. 
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The concept of transport driven by an external ratchet po- 
tential dates back to Pierre Curie J2]. He states that a locally 
asymmetric but globally symmetric system may induce global 
transport if it is kept out of equilibrium. Practical realizations 
include colloidal particles that may move either through peri- 
odic asymmetric micropores when driven by an imposed os- 
cillating pressure field J2l or they can be driven by a saw- 
tooth dielectric potential that is periodically switched on and 
off . Many different variations of ratchet mechanisms are 



nowadays studied I 
discrete objects 



. They are employed to transport or filter 
or to induce macroscopic transport of a 
continuous phase in systems without a macroscopic gradient 
(H |H 0] . Most modelling effort is focused on the former, but 
first models do as well exist for the latter |0, 0] . 



An experiment that at first sight seems unrelated has re- 
cently shown that drops may climb an inclined homogeneous 
substrate if it is vibrated harmoniously in a vertical direction 
IToll . The experiment is quite remarkable, as previously it had 
only been shown that substrate vibrations can 'unlock' drops 
pinned by substrate heterogeneities and therefore facilitate di- 
rected motion of drops in a global gradient. In particular, they 
allow to overcome effects of contact angle hysteresis on non- 
ideal substrates with a global wettability gradient fill [l2tl . A 
recent extension shows that drops can as well be driven by 
simultaneous vertical and horizontal substrate vibrations that 
are phase-shifted lfl3ll . Several hypothesis have been put for- 
ward as to why the vibrations induce the drop motion. Con- 
tact angle hysteresis, nonlinear friction, anharmonicity of the 
vibrations, convective momentum transport are all mentioned 
as possible 'minimal ingredients' fiol Il3l [Till . Anharmonic 
lateral vibrations with non-zero mean displacement are known 
to drive drop motion on horizontal substrates t\3\ . In general, 
several studies have adressed drops and free surface films on 
oscillating substrates since early work by Faraday, Kelvin and 
Rayleigh iflrjll . However, most consider fixed contact lines 



lfl7l [Till, purely inviscid flows lfl9l l20ll or a high frequency 
limit 12111 : and, most importantly, limit their study to vibra- 
tions either parallel or orthogonal to the substrate. 

In the present Letter we analyze a minimal hydrodynamic 
model for the situation depicted in Fig. Q] and show which 
of the above mentioned ingredients are not necessary for the 
drop motion to occur. We derive a thin film evolution equa- 
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FIG. 1: (color online) Sketch of a drop on a vibrating (frequency 
uj, vibration angle with the substrate normal /3) inclined substrate 
(inclination angle a). 

tion in the overdamped limit, i.e., based on a pure Stokes flow 
(no convective momentum transport) for a drop on an ide- 
ally smooth homogeneous inclined substrate (no contact an- 
gle hysteresis). The obtained results show that harmonic vi- 
brations are sufficient to drive a drop in a directed manner on 
a horizontal substrate lfL3ll or as well up an incline 

Oil- 
Analysing the underlying mechanism we find that the com- 
ponent of the harmonic vibration that is orthogonal to the sub- 
strate induces a nonlinear (anharmonic) response in the drop 
shape. As the latter determines the strongly nonlinear drop 
mobility this results in an asymmetric response of the drop to 
the vibration component that is parallel to the substrate. The 
induced symmetry breaking between forward and backward 
motion during the different phases of the vibration results in 
the observed net motion of the drop. This phenomenon can 
be seen as a rocked self -ratcheting of the drop Q4j] as it is the 
drop itself that introduces the local time-reflection asymme- 
try in the response to the time-periodic driving of the sliding 
motion. 

In the following we first present the model, then analyze it 
in the low frequency limit employing continuation techniques. 
Results for higher frequencies are obtained through numerical 
time integration. We use a long- wave approximation i22l l23ll 
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to describe the dynamics of a drop of liquid on an inclined 
substrate (inclination angle a) that is subject to harmonic vi- 
brations (frequency ui, at angle /3 to the substrate normal, see 
Fig.Q]). We consider a two-dimensional (2d) drop as we do not 
expect any conceptual difference to the 3d case. The resulting 
evolution equation for the film thickness profile h(x,t) is 



dth = d x 



377 



(d xP -f) 



(1) 



The divergence of the flow on the r.h.s. is expressed as the 
product of a mobility and the sum of a pressure gradient d x p 
and a lateral driving force /. -q is the dynamic viscosity. The 
pressure 



■ jd xx h - 11(h) + pgh [1 + a(t)] 



(2) 



contains the curvature pressure —jd xx h, where 7 denotes sur- 
face tension, the disjoining pressure H(h) that incorporates 
wettability [22, 24] and the hydrostatic pressure pgh[l+a(t)], 
where the time dependence results from the vibration compo- 
nent normal to the substrate. The lateral force 



/ = pg [a + pa(t)} 



(3) 



contains a constant part (force down the incline) and a time- 
modulated part (vibration component parallel to the sub- 
strate). The function a(t) — ao sin(wi) corresponds to the 
acceleration in units of g. The drop experiences a force of 
opposite sign, e.g., if the substrate is accelerated upwards 
and to the left the drop is flattened and pushed to the right. 
The partial wettability of the substrate is modelled through a 
precursor film model based on the disjoining pressure II = 
(A/h 3 + B/h 6 ) /6tt. It combines long-range destabilizing 
(A < 0) and short-range stabilizing (B > 0) van der Waals 
interactions lf25tl . 

Note, that consistency of the long-wave approach requires 
small free surface slopes as well as a, \f3\ <C 1 l23"tl . However, 
it is known that equations like ([TJ often predict the correct 
qualitative behavior even for systems with larger contact an- 
gles 1I22I1 . In consequence, we expect the here obtained results 
to hold qualitatively as well for larger (3. 

To nondimensional ize we introduce the scales to = 
Sjrj/hoK 2 , x — y/ 7/10 / and h — (B/A) 1 / 3 for t, 
x, and h, respectively, where n — A/Qnli^. The result- 
ing dimensionless equations correspond to Eqs. (JT]>-([3j with 
377 = 7 = B = 1 and A = — 1. The dimensionless modulated 
hydrostatic pressure and lateral force are ph = Gh [1 + a(t)] 
and / = G(a + f3a(t)) with G = pgho/n. The vibration 
period T = 2-k/uj is given in units of to- The fixed drop 
volume V = L(h — h p ) is determined by the domain size 
L, the mean film thickness h and the precursor film thick- 
ness hp = 1. The resulting transport along the substrate is 
measured after all transients have decayed and the vibration- 
induced shape changes of the drop are completely periodic in 
time. We quantify the transport by the mean drop velocity 
(v) = Ax/T where Ax denotes the distance the drop moves 
within one period T. 

Fig.|2]illustrates the typical behaviour of a drop during one 
vibration cycle. One notices forward/backward lateral motion 




FIG. 2: (color online) Space-time plot illustrating the evolution of 
the profile of a drop on a obliquely vibrated substrate during one 
vibration period. Each cycle results in the net motion of the drop to 
the left. The starting time is well after initial transients have decayed. 
Note that only part of the domain L is shown. Parameters are V = 
192, G = 0.001, (3 = 0.1, a = 0, T = 400, L = 128, a = 10. 



and small but significant changes of shape. In particular, the 
drop becomes flatter and broader during the first part of the 
cycle when its center of mass moves to the right (t < T/2, 
substrate acceleration is upward and to the left. In the second 
half of the cycle the drop becomes again higher and less wide 
while it moves to the left (t > T/2, substrate acceleration is 
downward and to the right). After one period the drop has 
moved a net distance to the left. In the course of one period, 
the orthogonal component of the oblique vibration modulates 
the hydrostatic pressure and causes a nonlinear response in 
the drop shape. That in turn determines the strongly nonlin- 
ear drop mobility [h 3 /3r/ in Eq. ([]])] and therefore induces 
an anharmonic response of the drop to the harmonic parallel 
vibration component that results in the observed net motion. 
Note, that a harmonic back and forth forcing alone results in 
zero net motion. 

First, we consider a slowly vibrating substrate, i.e. the typ- 
ical timescale of the intrinsic drop dynamics to is small com- 
pared to the vibration period T. In this low frequency limit the 
drop moves in a quasi-stationary manner, i.e., its shape and 
velocity at each instant during the vibration cycle are those 
of a stationary moving drop at the corresponding acceleration. 
The resulting family of moving drop solutions is parametrized 
by a(t) and can be obtained using continuation techniques 
ll2a 1271 1281. Averaging stationary drop velocities v(a(t)) over 
one vibration period gives the low frequency limit of (v). 

Figs. |3](a) and (b) show profiles of stationary moving drops 
and their velocities at various phases of the cycle on a horizon- 
tal substrate. As observed already in Fig. [2] during the first 
[second] half cycle the drop is compressed [decompressed] 
and slides to the right [left]. Fig.|4]shows the relationship be- 
tween acceleration a and the drop velocity v. The difference 
between the two respective curves for positive and negative 
acceleration is a measure of the anharmonic response of the 
drop. The mean velocity in the low frequency limit then cor- 
responds to twice the weighted area of the shaded region in 
Fig. |4{a) when taken from a = to a = ao- A numerical 



analysis gives (v) 



(v) - P and (v) ~ V 



1.67 



Note 
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FIG. 3: (color online) Shown are (a) several drop shapes and (b) the 
velocity during one vibration cycle for an obliquely vibrated hori- 
zontal substrate in the low frequency limit. For the used (3 = 0.3 the 
drop moves with v ~ —0.01, i.e., to the left. The blue dotted line 
in (b) indicates a harmonic variation of zero net flow. Remaining 
parameters are as in Fig. [2] 




FIG. 4: (color online) (a) Mean drop velocity depending on the accel- 
eration a for a horizontal substrate. The solid (dashed) line indicates 
a > 0, v > (a < 0, v < 0). (b) Phase diagram indicating the 
regions of sliding and climbing drops in the ao — /3~ 1/,2 -plane for a 
substrate inclination of a = 0.1. The results for (a) and (b) are in the 
low frequency limit with the remaining parameters as in Fig. [2] 



that for positive (3 in the present limit net transport is always 
directed to the left. 

From the transport properties on a horizontal substrate it is 
obvious, that such a vibration-caused motion can overcome a 
further external driving and, e.g., move a drop up an incline 
or against a wettability gradient, as long as the product /3a§ 
is larger than a critical value. This is confirmed by the phase 
diagram Fig.|4[b), as there the straight line ao ~ yVo/fi sepa- 
rates climbing from sliding drops, where vo is the drop veloc- 
ity for the given substrate inclination without any vibration. 
Note that the case j3 = a, i.e., a vibration vertical in the lab- 
oratory frame, is by no means special: uphill motion is gen- 
erated above a finite threshold acceleration. However, in the 
limit /3 — ► 0, i.e., for a vibration normal to the substrate, the 
threshold acceleration becomes infinite and the net transport 
is zero. 

For practical purposes, e.g., for microfluidic applications, 
drops need to be transported in a limited time, i.e., the be- 
haviour at sufficiently large frequencies has to be understood. 



FIG. 5: (color online) Master curve for the scaled mean velocity 
(v) I f3a,Q as a function of the vibration period T for sets (ao, /3) as 
given in the legend. The horizontal dashed line indicates the result in 
the low frequency limit. Remaining parameters are as in Fig. [2] 
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FIG. 6: (color online) Dependence of transport on drop volume: 



(a) shows the special vibration periods T n: 



V 1 



and 



T c ~ V , as detailed in the text; (b) gives the absolute value of 
the mean velocities (y) max ~ 

^2.78*0.06 at period ^ ^ an( j m [he 

low freqency regime (Too)- Remaining parameters are as in Fig. [2] 



We study that regime employing an adaptive time-step 4th or- 
der Runge-Kutta method and prevent numerical instabilities 
by a switching upwind differencing for the driving term. The 
above deduced scaling (v) ~ a^[3 still holds in a large fre- 
quency range, illustrated in Fig.|5]by collapsing dependencies 
of mean velocity on vibration period for various sets (ao, (3) 
onto a single master curve. A second interesting feature is 
the observed reversal of net transport at a small critical pe- 
riod T c ?s 55 that results in a (small) positive mean velocity 
even for a horizontal substrate (a = as in Fig. [5]). The re- 
versed flow is strongest at T max « T c /2 r* 27, but always 
about one order of magnitude smaller than the transport to the 
left in the low frequency regime. Note that in agreement with 
the first observation T c and T max are nearly independent of ao 
and (3. They do, however, depend on drop volume (see Fig.|6| 
indicating that the flow reversal might be triggered when the 
vibration frequency becomes larger than an eigen frequency 
of the drop, i.e., for larger imposed frequencies the response 
of the drop becomes delayed and phase-shifted w.r.t. the forc- 
ing. This is corrobated by an inspection of drop profiles over 
the course of a cycle (not shown). This result is similar to the 
dependence of the direction of net motion on the phase-shift 
between vertical and horizontal vibration observed in lfl3ll . 
Finally, Fig. [7] presents a phase diagram that shows the 
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FIG. 7: (color online) Phase diagram of sliding and climbing drops 
in the (an, l/T) plane for an obliquely vibrated inclined substrate. 
Beside a = 0.05 parameters are as in Fig. [2] 

influence of the vibration period and peak acceleration on 
the transport direction for fixed substrate inclination angle 
a = 0.05. The critical acceleration a c (T _1 ) separating slid- 
ing and climbing drops approaches a constant at small fre- 
quencies and diverges at finite T~ . As we do not model 
any effects of non-ideal substrates (contact angle hysteresis, 
contact line pinning) our phase diagram differs in two aspects 
from the experimental one ifioll : (i) we find a monotonic de- 
crease of a c with decreasing frequency (T" 1 ) whereas ifioll 
finds evidence for a small increase at small T -1 ; and (ii) 
here drops always slide down below a c whereas in the experi- 
ments one finds instead a transition between static and climb- 
ing drops at large T -1 . 

To conclude, we have proposed and analyzed a minimal hy- 
drodynamic model for the experimentally observed drop mo- 
tion that is driven by harmonic oblique substrate vibrations 



ifioll . Our analysis has ruled out convective momentum trans- 
port, anharmonicity of vibrations and contact angle hysteresis 
as necessary for the motion. The mechanism that moves drops 
on a horizontal substrate or even up an incline is based on a 
nonlinear response of the drop shape to the normal vibration 
component. This breaks the back-forth symmetry of the re- 
sponse of the drop to lateral oscillations and therefore causes 
a net motion. The found mechanism is in line with the hy- 
pothesis put forward in ifioll based on a mechanical analogue 
that a breaking of the front-back symmetry due to a nonlinear 
friction law is sufficient to induce transport. Here, we have 
identified the strongly shape-dependent nonlinear mobility in 
Eq. (Q]i as the relevant hydrodynamic 'nonlinear friction'. Be- 
side establishing the basic mechanism, our analysis has re- 
vealed several interesting features that should be investigated 
in future experiments. In particular, we have identified a num- 
ber of scaling laws valid in a broad frequency range and found 
a flow reversal at high frequencies. 

We have argued that the phenomenon can be seen as a self- 
ratcheting of the drop as its shape changes are instrumental 
in producing local time-reflection asymmetries. Comparing 
Eq. ([T]i with Fokker-Planck descriptions for ratchet systems 
of interacting particles @ |29| one may further the analogy by 
noting that (i) the lateral vibration component corresponds to 
an imposed rocking, (ii) the orthogonal vibration component 
corresponds to an imposed temporal temperature modulation, 
(iii) the role of the spatial asymmetry in the ratchet potential is 
here taken by non-linear couplings due to strongly nonlinear 
prefactors of the diffusion (2nd order) and the transport (1st 
order), (iv) our surface tension term is analogous to a mean 
field expansion of the distribution function up to second order. 
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